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Quantum transport for DNA conduction has widely studied with interest in application as a 
candidate in making nanowires as well as interest in the scientific mechanism. In this paper, we 
review recent works with concerning the electronic states and the conduction/transfer in DNA 
polymers. We have mainly investigated the energy band structure and the correlation effects of 
localization property in the two- and three-chain systems (ladder model) with long-range correlation 
as a simple model for electronic property in a double strand of DNA by using the tight-binding model. 
In addition, we investigated the localization properties of electronic states in several actual DNA 
sequences such as bacteriophages of Escherichia coli, human-chromosome 22, compared with those 
of the artificial disordered sequences with correlation. The charge transfer properties for poly(dA)- 
poly(dT) and poly(dG)-poly(dC) DNA polymers are also presented in terms of localization lengths 
within the frameworks of the polaron models due to the coupling between the charge carriers and 
the lattice vibrations of the double strand of DNA. 

PACS numbers: 



I. INTRODUCTION 

Recent interests on semiconducting DNA polymers 
have been stimulated by successful demonstrations of the 
nanoscale fabrication of DNA, where current- voltage (I- 
V) measurements for poly(dA)-poly(dT) and poly(dG)- 
poly(dC) DNA polymers have been done [1it4|. For such 
artificial periodic DNA systems, the energy band struc- 
ture is a useful starting point in order to interpret the 
experimental results such as semiconductivity and the 
metal-insulator transition 

On the other hand, Tran et al. measured conductivity 
along the double helix of lambda phage DNA (A— DNA) 
at microwave frequencies, using the lyophilized DNA in 
and also without a buffer [6] . The conductivity is strongly 
temperature dependent around room temperature with a 
crossover to a weakly temperature dependent conductiv- 
ity at low temperatures [3] . Yu and Song showed that the 
A— DNA can be consistently modeled by considering that 
electrons may hop through the variable range hopping for 
conduction without invoking additional ionic conduction 
mechanisms, and that electron localization is enhanced 
by strong thermal structural fluctuations in DNA 7]. 
Indeed, the sequence of base pairs(bp) of the A— DNA 
is inhomogeneous as in disordered material systems. 

Moreover, charge migration in DNA have been mainly 
addressed in order to clarify the mechanism of damage 
repair which are essential to maintain the integrity of 
the molecule. The precise understanding of the DNA- 
mediated charge migration would be important on the 
descriptions of damage recognition process and protein 
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binding, or on the engineering biological processes [8|. 
The stacked array of DNA bp provides an extended path 
to a long-range charge transfer although the dynami- 
cal motions or the energetic sequence- dependent hetero- 
geneities are expected to reduce the long-range migra- 
tion. Photoexcitation experiments have unveiled that 
charge excitations can be transfered between metalloin- 
tercalators through the guanine highest occupied molec- 
ular orbitals of the DNA bridge t^]. The subsequently 
low-temperature experiments showed that the radiation- 
induced conductivity is related to the mobile charge car- 
riers, migrating within frozen water layers surrounding 
the DNA helix, rather than through the base-pair core. 

Those experiments are summarized as follows: (i) 
Band gap reduction of a double strand of DNA, (ii) 
Transition from the tunneling hopping to the band hop- 
ping, (iii) Anomalously strong temperature dependence 
of band gap, (iv) Highly nonlinear temperature depen- 
dence of the DC conductivity, (v) Low conductivity of 
DNA with a complicated sequence such as A— DNA and 
high conductivity of DNA with a simple sequence such 
as poly(dG)-poly(dC) and poly(dA)-poly(dT). These re- 
sults suggest that the anomalously strong temperature 
dependence on the physical quantities is attributed to 
the "self-organised" extrinsic superconductive character 
of DNA due to the formation of donors and accepters. 
Here we would like to note the following: Usually, as 
in solid state physics, the extrinsic semiconductor is re- 
alized when external impurities are introduced in pure 
substrate materials. Such impurities produce the extrin- 
sic semiconductive nature. However, in real DNA helices, 
there are no such impurities from outside; but there ex- 
ist already a complicated arrangement of bases of adenine 
(A), guanine(G), cytosine (C) and thymine (T) inside of 
the DNA, where among the bases each one of bases may 
regard other bases as impurities. Hence, a kind of "self- 
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organized" extrinsic semiconductive nature may appear. 

However, the electronic transport properties in DNA 
are still controversial mainly due to the complexity of the 
experimental environment and the molecule itself. Al- 
though theoretical explanations for the phenomena have 
been tried by some standard pictures used in solid state 
physics such as polarons, solitons, electrons and holes, 
the situation is still far from unifying the theoretical 
scheme. 

Each DNA sequence is packed in a chromosome, vary- 
ing in length from lO^bp in yeast to lO^bp in human. 
In general, the length of a mutation is relatively short 
(106p) as compared to the length of a gene (10'^ — lO^bp). 
Because the mutation rate is very law the mechanical 
and thermodynamic characteristics are maintained for 
the mutation |20|]. In Sect. 6, we give a brief discussion 
about the mutation as the proton transfer between the 
normal and tautomer states. 

The Watson-Crick(W-C) base-pair sequence is essen- 
tial for DNA to fulfill its function as a carrier of the 
genetic code. The specific characteristics is also used 
extensively even in various fields such as anthropology 
and criminal probe. As observed in the power spectrum, 
the mutual information analysis and the Zipf analysis 
of the DNA base sequences such as human chromosome 
22(IICh-22), the long-range correlation exists in the to- 
tal sequence, as well as the short-range periodicity. In 
this article, we mainly discuss the relationship between 
the correlation in the DNA base sequences and the elec- 
tronic transport /transfer. 

The electrical transport of DNA is closely related 
to the density of itinerant tt— electrons because of the 
strong electron-lattice interaction. Resistivities of two 
typical DNA molecules, such as poly(dG)-poly(dC) and 
A— DNA, are calculated. At the half-filling state, the 
Peierls phase transition takes place and the poly(dG)- 
poly(dC) and poly(dA)-poly(dT) DNA polymers ex- 
hibit a large resistivity. When the density of itinerant 
TT— electrons departs far from the half-filling state, the 
resistivity of poly(dG)-poly(dC) becomes small. For the 
A— DNA, there is no Peierls phase transition due to the 
aperiodicity of its base pair arrangement. The resistiv- 
ity of poly(dG)-poly(dC) decreases as the length of the 
molecular chain is increasing, while that of A— DNA in- 
creases as the length is increasing. 

In Sect. 3, we introduce the Hiickel model of DNA 
molecules to treat with the tt— electrons. Moreover, in 
Sect. 4, we give the effective polaron model including the 
electron-lattice coupling dynamics. In Sect. 5, we present 
correlation effects of localization in the ladder models and 
the formation of localized polarons (Holstein's polarons 
or solitons) due to the coupling between the charge car- 
riers and the lattice vibrations of the double strand of 
DNA. These polarons in DNA act as donors and accep- 
tors and exhibit an extrinsic semiconductor character of 
DNA. The results are discussed in the context of experi- 
mental observations. 

The present review article essentially follows the line 



in Refs. [TT| - [l5| . In particular, it is written with a 
view of the localization and/or delocalization problem 
in the quasi-one-dimensional tight-binding models with 
disorder. In appendices we give some calculations and 
explanations related to main text. 



II. CORRELATED DNA SEQUENCES 

As is well known, DNA has the specific binding proper- 
ties, i.e., only A-T and G -C pairs are possible, where the 
bases of nucleotide are A, T, G, C. The backbones of the 
bases, sugar and phosphate groups, ensure the mechani- 
cal stability of the double helix and protect the base pairs. 
Since the phosphate groups are negatively charged, the 
topology of the duplex is conserved only if it is immersed 
into an aqueous solution containing counterions such as 
Na'^ and Afg"*" that neutralize the phosphate groups. 

The clustering of similar nucleotides can be clarified by 
studying the properties of the cluster size distributions 
on the various real DNA sequences, ranging from the 
viral to higher eukaryotic sequences. It is shown that 
the distribution function P{S) about the number S of 
the consecutive C-G or A-T clusters becomes P{S) oc 
exp{— aS*} [H, The values of the scaling exponent a 
of CG are much larger than a of AT. The maximum value 
of the A-T cluster size is found to be much larger than 
that of the C-G cluster size, which implies the existence 
of large A-T clusters. 

Moreover, it has been found that the base sequences of 
various genes exhibit a long-range correlation, character- 
ized by the power spectrum S{f) ~ (0.1 < /3 < 0.8) 
in the low frequency limit (/ << 1) [I^. As was observed 
in the power spectrum, the mutual information analysis 
and the Zipf analysis of the DNA base sequences such as 
the HCh-22, the long-range structural correlation exists 
in the total sequence as well as the short-range period- 
icity. The eukaryote's DNA sequence has an apparently 
periodic repetition in terms of the gene duplication. The 
correlation length in the base sequence of genes changes 
from the early eukaryote to the late eukaryote as a result 
of evolutionary process. It is found that the long-range 
correlation tends to manifest in the power spectrum of 
the total sequence rather than in the power spectra of 
the exon and intron parts, separately [19|- 

On the other hand, in a DNA molecule, the charge 
carriers move along a double-helix formed by two com- 
plementary sequences of four basic nucleotides: A, T, 
G, and C. A conduction band would form, if the DNA 
texts would exhibit some periodicity. The electrical re- 
sistance of the DNA molecule strongly fiuctuates even if 
a single nucleotide in a long sequence is replaced (or re- 
moved). Quantum transport through the DNA molecule 
is also strongly affected by the correlations. The localiza- 
tion property of a single-chain disordered systems with 
long-range correlation has been also extensively studied. 
The correlated disorder can lead to delocalized states 
within some special energy windows in the thermody- 
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namic limit. 

Accordingly, it is very interesting to compare the lo- 
calization nature of the electronic states in the real DNA 
sequence with those in the artificial disordered sequence 
with long-range correlation. Recently, Krokhin et al. 
have reported that much longer localization length has 
been observed in the exon regions than in the intron 
regions for practically all the allowed energies and for 
all randomly selected DNA sequences [l^]. Through the 
statistical correlations in the nucleotide sequence, they 
suggest that the persistent difference of the localization 
property is related to the qualitatively different informa- 
tions stored by exons and introns. 

In the Sect. 5, we numerically give localization nature of 
the electronic states in some real DNA sequences such as 
bacteriophages of Escherichia coli (E.coli) and HCh-22, 
and so on in ladder models HH . We also investigate 
the correlation effect on the localization property of the 
one-electronic states in the disordered ladder models with 
a long-range structural correlation that is generated by 
the modified Bernoulli map. Obviously, the correlation in 
the DNA sequences affects not only the electronic con- 
duction but also the twist vibration of the backbones. 
However, such effects can be approximately ignored in 
the scope of this article. 



III. HOMO-LUMO GAPS 

Each compositional unit of a DNA polymer is complex 
although the DNA polymers can be regarded as quasi- 
one-dimensional systems. (See Fig. 1.) In this section, we 
give a brief review on the Hitckel theory in order to ana- 
lyze the relationship between electronic structures of the 
separate nucleotide groups and of their infinite periodic 
chains. Ladik tentatively concluded that electrons hop 
mainly between the bases along the helical axis of DNA, 
so that it is enough to take into account the overlap inte- 
grals between the tt— orbitals of the adjacent base pairs in 
the DNA duplexes [1]. These studies have revealed that 
DNA polymers are insulators with an extremely large 
band gap Eg (about 10-16 eV) and narrow widths of the 
valence and conduction bands (about 0.3-0.8 eV). This is 
a consequence of the orbital mixing between the highest 
occupied molecular orbital(HOMO) and the lowest unoc- 
cupied molecular orbital (LUMO) in the base groups of 
nucleotide in the DNA. We can control the charge injec- 
tion into DNA and the conduction in DNA, by adjusting 
the HOMO-LUMO of the polymer and the Fermi energy 
of electric leads, respectively. 

Brumel et al. found that semiconduction in DNA is 
very important for such systems since the activation en- 
ergies of nucleotides are lower than those of nucleosides 
[2^ . This was the first suggestion that propagation along 
the sugar and phosphate groups play a significant role in 
the transport properties of DNA as well as propagation 
along the base-stacking of the nucleotide groups in the 
center of the DNA molecule. 




FIG. 1: The single bases of A, G, C, T. Here dots (•) mean 
TT-electrons. 



We revisit the problem of the electronic properties of 
individual molecules of DNA, in order to know electron 
transport in the double strand of DNA as a mother ma- 
terial for the single and double strand of DNA, by taking 
into account only tt— electrons in the system 5]. To do 
so, we review the theory of tt— electrons in DNA, using 
the Hitckel approximation for tt— electrons in both the 
sugar-phosphate backbone chain and the tt— stacking of 
the nitrogenous bases of nucleotide. 

A. Huckel Model 

Applying the basic knowledge of quantum chemistry to 
biomolecules of nitrogenous bases, we can find the total 
number of tt— orbitals and tt— electrons in the bases. It is 
summarized in Table 1. 

Let us consider the famous Hiickel model in quan- 
tum chemistry (22| - [27| . This model concerns only the 
TT— orbitals in the system. In this context, this theory 
is closely related to the so-called tight-binding model 
in solid state physics, which concerns very localized or- 
bitals at atomic sites such as the Wannier's wavefunction. 
Therefore, this approach has been extensively applied to 
many polymer systems such as polyacetylene with a great 
success. In the so-called Hitckel approximation we adopt 
the orthogonality condition for the overlap integrals 

Srr = 1, Sr, = [r ^ s), (1) 

where Srs denotes overlap integrals between atomic or- 
bitals at rth and sth sites. We assume the special form 
of the resonance integrals 

hrs — S'rs{hf"r -\- /iss), (2) 
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TABLE L The total numbers \N\, \C\, \0\ and |P| of N, C, O 
and P atoms and of 7r-orbitals and 7r-electrons in the bases of 
A, G, C and T, and the sugar(S) and phosphate(Ph) groups, 
respectively. Here there are 12, 14, 10 and 10 7r-electrons for 
the 10, 11, 8, 8 TT-orbitals in the A, G, C and T base molecules, 
respectively, while there are 8 and 8 7r-electrons for the 4 and 
5 TT-orbitals in the sugar and phosphate groups, respectively. 
Here we note that the carbon atom of CH3 in the T molecule 
does not have any 7r-electron since it forms the sp'^— hybrid 
orbitals. This provides 8 7r-electrons for the T molecule. 

with K — 1.75 and the overlap integrals Srs {r = s) are 
not necessary to be diagonal, and otherwise. Usually the 
on-site (r = s) resonance integrals are called the Coulomb 
integrals denoted by a^, while the off-site (r 7^ s) inte- 
grals are called the resonance integrals denoted by f3rs 
such that ar — hrr, Prs — hrs {r 7^ s). Here we would 
like to emphasize the following. In the sense of the Hitckel 
theory, the parameters are taken empirically. This means 
that the parameters are adjustable and feasible to give 
consistent results with the experimental results or the ab 
initio calculation results. Therefore, the exact values of 
the parameters are neither so important nor should be 
taken so seriously in this framework. This is because 
once one can obtain much more precise values for the 
Hitckel parameters, one can provide the more plausible 
results from the Hiickel theory. Although many efforts 
of the ab initio calculations have been done for DNA 
systems, unfortunately at this moment there seem to be 
very few first principle calculations for such parameters 
in the DNA systems to fill out this gap. Nevertheless, 
we must assign some values for the H-iickel parameters 
in order to calculate the electronic properties of DNA in 
the framework of the Hiickel theory. So, we look back to 
the original method about time when the Hiickel theory 
was invented. 

For this purpose to use the standard Hiickel theory, let 
us adopt some simple formulae for the Hiickel parame- 
ters, which are defined as follows: Let X and Y be two 



different atoms. Denote by ax the Coulomb integral at 
the X atom and by I^xy the resonance integral between 
the X and Y atoms: 

ax = a + ax (3, Pxy = IxyP- (3) 

Here ax and txY are the empirical parameters that are 
supposed to be adopted from experimental data. And the 
parameters a and /? are important. These can be thought 
of the fundamental parameters in our problem of biopoly- 
mers. Conventionally we take a as the Coulomb integral 
for the 2^2;— orbital of carbon and /3 as the resonance 
integral between the 2pa:— orbitals of carbon, such that 
a = ac = 0, (3 — (3cc = 1- This means that the energy 
level of a carbon atom is taken as the zero level, and the 
energy is measured in units of the resonance integral be- 
tween carbon atoms. We note that the empirical values 
obtained from experiments are usually given by 

a w -6.30 - -6.61ey, /3 w -2.93 - -2.95el/. (4) 

Since the biomolecules consist of the atoms C, N, O, P, let 
us find the plausible values of the Hiickel parameters for 
them to apply the Hiickel model to biomolecules of DNA. 
The Pauling's electronegativity for carbon (C), nitrogen 
(N), oxygen (O) and phosphorus (P) are the following: 

Xc ~ 2.55, XN « 3.0, xo ~ 3.5, xp ~ 2.1. (5) 

Then we can defined the Hiickel parameters for carbon, 
using the Sandorfy's form ula [28| . Mulliken's formula[29j 
and Streitwieser's formula[30|, we find 
a + 0.7/3, = a + 1.53/3, ap = a- 0.724/3. and we ob- 
tain = a -I- 2.53/3, a^ ~ a 0.276/3. and /3p_^ = /3, 
Pc^c = 1-1/3. Pc^N = Pc-N = 0.8/3, /3c^^ = 1.1/3, 
fi(j_Q ~ 0.9/3, /3(^^Q = 1.7/3. Here a^(aj^) denotes the 
Coulomb integral of atomic state X with one(two) elec- 
tron(s) occupied. See appendix A for the formula. And 
if one can get more accurate values from the ab initio 
calculations, then we can always replace the Coulomb 
integrals by the new set of values. 



B. HOMO-LUMO of Biomolecules 

For example, considering the topology of the hopping 
of electrons on the 7r-orbitals, we now find the following 
Hiickel matrices Ha for A: 
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Pauli's exclusion principle tells us that each state with an 
energy level is occupied by a pair of electrons with spin 
up and down. So, tt— electrons occupy the energy levels in 
the spectrum from the bottom at low temperature. Since 
the lower energy levels with one half of the total number 
of TT— electrons can be occupied by the tt— electrons, there 
appears an energy separation between the occupied and 
the unoccupied states, which is called the energy gap. 

The energy levels of the HOMOs and LUMOs are given 
by £h = 0.888, El — —0.789 for A, where the energy is 
measured in units of /?. We can do the same calcula- 
tion for G,C,T, respectively. Defining the energy gap 
between the LUMO and HOMO, Ae = - e^- Then 
we obtain the result as, Ae^i = 1.677, Aeg — 1.555, 
Aec = 1.535, Aet = 1-713. (See Fig.2(a).) The to- 
tal energy Etot = 2 Y.j=occ states ^3 7r-electrons 
of A, G, C, T are Etot{A) = 20.74, Etot{G) = 26.47, 
Etot{C) = 19.05, Etot{T) = 21.14, respectively. There- 
fore, 

Etot{C) > Etot{A) > EtotiT) > Etot{G). 

This shows that since the lower the ground state energy 
the more stable the system, the most stable molecule is 
G while the most unstable molecule is C. 

HOMOs and LUMOs of sugar and phosphate may 
be required when we consider charge conduction of the 
DNA polymer because the sugar and phosphate group 
constituting the backbones of DNA polymer also have 
TT— electrons. We give the results for the single sugar- 
phosphate in Fig. 2(b). Moreover, the results based on 
the Hiickel approximation have been extended to the sin- 
gle nucleotide systems such as A, G, C, T with the single 
sugar-phosphate group, and the system of a single strand 
of DNA with an infinite repetition of a nucleotide group 
such as A, G, G and T, respectively. See Ref.([5j) for 
the more details. This reorganization, which is difficult 
to calculate due to the complexity of the combined sys- 
tem, may lead to a smaller HOMO-LUMO gap and wider 
band widths than for the bare molecule. 

When the system of the nucleotide bases such as A, G, 
G, T exists as an individual molecule, there is always an 
energy gap between the LUMO and HOMO states, where 
the order of the gap is about several eV. This means that 



(a) : 




FIG. 2: (a) The spectrum of the 7r-orbitals of A, G, C, T. (b) 
The spectrum of the 7r-orbitals of a sugar-phosphate group. 
PO4 stands for the phosphate group where the electron hop- 
ping between the oxygen sites is taken into account, while 
PO'4 stands for the the phosphate group where the electron 
hopping between the oxygen sites is forbidden. SP stands for 
the sugar-phosphate group where the electron hopping be- 
tween the oxygen sites is taken into account. The energies 
are measured in units of j/3|. Here dots (•) mean 7r-electrons 
and the level with four dots means the double degeneracy of 
the level. 



the nucleotides of A, G, C, T have the semiconducting 
character in its nature. If the tt— stacking of the base 
is perfect, then there are two channels for tt— electron 
hopping: One is the channel through the base stacking 
and the other is the one through the backbone chain of 
the sugar-phosphate. In this limit, the localized states of 
the original nucleotide bases become extended such that 
the levels form the energy bands. 

Thus, we believe that our Hiickel approach in this pa- 
per fills out the gap between the simple approach of 
mathematical models for tight-binding calculations and 
the approach of the quantum chemical models for ab ini- 
tio calculations from the first principle, where in the for- 
mer we assume one orbital with one electron at one nu- 
cleotide while in the latter we include all electrons and 
atoms in the system. 
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IV. EFFECTIVE POLARON MODELS 

First principle methods are powerful enough to under- 
stand the basic electronic states of the molecules. On 
the other hand, the complementary model-based Hamil- 
tonian approach is effective for understanding those of 
polymers as well. In this section, we give some effective 
ID models and the basic properties of polymers such as 
trans-polyacetylene and DNA polymers. In the last sec- 
tion, we treated HOMOs and LUMOs and the occupation 
of orbitals by the tt— electrons with spin. Hereafter, we 
omit spin of electrons when we consider the one-electron 
problems, for the sake of simplicity in notation. 

A. Tight-binding Models for Polymers 

It is known that the electronic properties of planar 
conjugated systems are dominated by tt— systems with 
one orbital per site. The electronic Hamiltonian for 2pz 
orbitals through a tight-binding model with the nearest 
neighbors interactions only: 

— EnC\Cn — Kiri+l(C'iCn+l + C'„+iC^OQ) 

n n 

where Cn and Cj^ are creation and annihilation operators 
of an electron at the site n. The matrix elements are 



where C^{C„) creation (annihilation) operator of the car- 
rier at the site n 0, [U, HI]. The carrier site energies 
En are chosen according to the ionization potentials of 
respective DNA bases as, ea — 8.24eV, er — Q.lAeV, 
ec = 8.87ey, and eq — 7.75eV, while the hopping in- 
tegral, simulating the tt — tt stacking between adjacent 
nucleotides, is taken as tuNA = 0.1 — OAeV. 9n,n+i 
denotes the relative twist angle deviated from equilib- 
rium position between the nth and (n -I- l)th molecules 
due to temperature T. Then we can estimate the or- 
der of the hopping integral. Let 9n,n+i is an indepen- 
dent random variable that follows Gaussian distribution 
with {0n,n+i) = and (On^n+i) = kBT/Iujj, where / is 
the reduced moment of inertia for the relative rotation 
of the two adjacent bases and w/ is the oscillator fre- 
quency of the mode {luij = 250K) . Therefore the fluctu- 
ation of the Hopping term is about tDNA{cos{9n,n+i)) ~ 
toNAi^ — (^?i.n+i)/2) ~ O.lGeV at room temperature, 
which is much smaller than that lAeV) in the diag- 
onal part. Accordingly, as a simple approximation, we 
can deal with the diagonal fluctuation with keeping the 



obtained from the extended Hiickel theory as given in 
the previous section: 



En = -a„ (7) 
Km = (A72)(a„ -Ka™)5„„, (8) 



where a„ is ionization energy (Coulomb integral) of the 
nth orbital, and Snm is overlap integral (resonance 
integral) between the nth and mth orbitals centered at 
neighbour sites given in Sect. 3. We usually treat the 
copolymer as a tt— system with one orbital per site and 
represent the electronic Hamiltonian for the 2pz carbon 
and nitrogen orbitals through a tight-binding model with 
nearest neighbor interactions. DNA polymers are also 
believed to form an effectively one-dimensional molecular 
wire, which is highly promising for diverse applications. 
Basically, the carriers mainly propagate along the aro- 
matic TT — TT stacking of the strands (the intcrstrand cou- 
pling being much smaller), so that the one-dimensional 
tight-binding chain model can be a good starting point 
to minimally describe a DNA wire[? ]. 

Bruinsma et al. also introduced the effective tight- 
binding model that describes the site energies of a carrier 
located at the nth molecule as. 



(9) 

I 

hopping integral constant. 



Roche investigated such a model for poly(GC) DNA 
polymer and A— DNA, with the on-site disorder arising 
from the differences in ionization potentials of the base 
pairs, and with the bond disorder tjjNA cos(6n,n+iit)) re- 
lated to the random twisting fluctuations of the nearest- 
neighbor bases along the strand [s^. While for poly(GC) 
the effect of disorder does not appear to be very dramatic, 
the situation changes when considering A— DNA. 



Actually, modern development of physico-chemical ex- 
perimental techniques enables us to measure directly 
the DNA electrical transport phenomena even in single 
molecules. Moreover, several groups have recently per- 
formed numerical investigations of localization properties 
of the DNA electronic states based on the realistic DNA 
sequences. 
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n n 
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B. Twist Polaron Models for DNA Polymers 

In this subsection we give an effective ID model with 
realistic parameters for the electronic conduction of the 
poly(dA)-poly(dT) and poly(dG)-poly(dC) DNA poly- 
mers. 

Specially, a distinctive feature of biological polymers 
is a complicated composition of their elementary sub- 
units, and an apparent ability of their structures to sup- 
port long-living nonlinear excitations. In their polaron- 
like model, Hennig and coworkers studied the electron 
breather propagation along DNA homopolynucleotide 
duplexes, i.e. in both poly(dG)-poly(dC) and poly(dA)- 
poly(dT) DNA polymers, and, for this purpose, they esti- 
mated the electron- vibration coupling stren gth in DNA, 
using semiempirical quantum chemistrv |33l435l | . Chang 
et al. have also considered a possible mechanism to ex- 
plain the phenomena of DNA charge transfer. The charge 
coupling with DNA structural deformations can create a 
polaron and thus promote a localized state. As a re- 
sult, the moving electron breather may contribute to the 
highly efficient long-range conductivity [s^. Recent ex- 
periments seem to support the polaron mechanism for 
the electronic transport in DNA polymers. 

The Hamiltonian for the electronic part in the DNA 
model is given by 

H,,{t) = En{t)ClCn - K„+i(i)(C,tC„+i + C„+ 

n n 

where C„ and Cj^ are creation and annihilation operators 
of an electron at the site n. The on-site energies En{t) 
are represented by 

E,,{t) = Eo + kr„{t), (11) 

where Eq is a constant and r„ denotes the structural 
fluctuation caused by the coupling with the transver- 
sal Watson-Crick H-bonding stretching vibrations. The 
schematic illustration is given in Fig. 3. The transfer in- 
tegral Vnn+i depends on the three-dimensional distance 
dnn+i between adjacent stacked base pairs, labeled by n 
and n -I- 1, along each strand. And it is expressed as 

Vnn+l{t) = Fo(l-ad„„+l(t)). (12) 

Parameters k and a describe the strengths of interac- 
tion between the electronic and vibrational variables, re- 
spectively. The 3D displacements dnn+i also give rise 
to variation of the distances between the neighboring 
bases along each strand. The first order Taylor expansion 
around the equilibrium positions is given by 

dnn+l{t) = -^(l-cos0o)(r„(t)-Hr„+i(t)). (13) 

i?o represents the equilibrium radius of the helix, 9q is the 
equilibrium double-helical twist angle between base pairs, 
and ^o the equilibrium distance between bases along one 
strand given by 

4 = [a' +mlsxn\e^/2)Y'\ (14) 
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kAT 


0.778917 eVA-' 


CtAT 


0.053835 A"^ 


kcc 


-0.090325 eVk-^ 




0.383333 A"^ 



TABLE II: Basic parameters for DNA molecules. The 
subscripts, AT and GC, for k and a denote for ones of 
the poIy(dA)-poly(dT) and poIy(dG)-poly(dC) DNA poly- 
mers, respectively. The other parameters are _Eo = Q.leV , 
Vo = O.leV, a = 3.4A, i?o = lOA and Bo = 36°. 



with a being the distance between the neighboring base 
pairs in the direction of the helix axis. We adopt real- 
istic values of the parameters obtained from the semi- 
empirical quantum-chemical calculations. (See table 2.) 
Further, we consider {?"„} as independent random vari- 
ables generated by a uniform distribution with the width 
(r„ e [— M^, VF]). Accordingly, fluctuations in both the 
on-site energies and the off-diagonal parts in the Hamil- 
tonian (1) are mutually correlated because they are gen- 
erated by the same random sequence r„. (See Fig. 2(c).) 
The typical value of ly is ly = 0.1[A], which approxi- 
mately corresponds to the variance in the hydrogen bond 
lengths in the Watson-Crick base pairs, as seen in the X- 
^ll^y diffraction experiments [13 . 

A quasi-continuum spectrum is possessed of a wealth 
of dynamical modes. In principle, each of these can in- 
fluence the DNA charge transfer/transport. But, since 
there are more or less active modes, it is possible to 
take the whole manifold of the DNA molecular vibra- 
tions into two parts: vibrations which are most active 
and a stochastic bath consisting of all the other ones. 
Computer simulations have pinpointed that the dynam- 
ical disorder is crucially significant for the DNA trans- 
fer/transport. Several attempts to formulate stochasti- 
cal models for the interplay of the former and the lat- 
ter have already appeared in the literature(see, for ex- 
amples). In this paper, we will deal with the polaron- 
like model of Hennig and coworkers as described in the 
works [s^. [35} ■ where charge-|-breather propagation along 
DNA homopolynucleotide duplexes [i.e. in the both 
poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA poly- 
mers] has been studied. 



C. Localization due to Static Disorder 

In this subsection, we assume that the fluctuations of 
{r„} are frozen (quenched disorder) and independent for 
different links; that is, we investigate localization prop- 
erties of the poly(dA)-poly(dT) and poly(dG)-poly(dC) 
DNA polymers 3q. 

The Schrodinger equation H^i]^) = is written in 



8 




FIG. 3: Sketch of the structure of the DNA model. The bases 
are represented by bullets and the geometrical parameters 
Ro,io,do,rn+i and d„„+i are indicated. 



the transfer matrix form, 



E-E„ 



1 







(15) 



where 0„ is the amphtude of the electronic wavefunction 
1$) = X]ra^"l")> where \n) — Cj^|0) at the base pair 
site n. |0) is the Fermi vacuum. We use the localization 
length ^ and/or Lyapunov exponent 7 calculated by the 
mapping ()15|) in order to characterize the exponential lo- 
calization of the wave function. Originally the Lyapunov 
exponent (the inverse localization length) is defined in 
the thermodynamic limit (iV — > 00). However, here we 
use the following definition of the Lyapunov exponent for 
the electronic wave function with a large system size N 

MM. 



'N- 



2N 



(16) 



We use the appropriate initial conditions (/)o = 0i = 1, 
and for large 7V(>> ^) the localization length and the 
Lyapunov exponent are independent of the boundary 
condition. The energy- dependent transmission coeffi- 
cient T{E, N) of the system between metallic electrodes 
is given as T{E, N) = exp{—2'yN) and is related to Lan- 
dauer resistance via g = (1 — r)/Tin units of the quan- 
tum resistance /i/2e2(~ 13[fcfim]) [43,|4l|. 

We have numerically investigated the localization 
properties of electronic states in the adiabatic polaron 
model of poly(dG)-poly(dC) and poly(dA)-poly(dT) 
DNA polymers with the realistic parameters obtained us- 
ing the semi-empirical quantum-chemical calculations. 

We compare the localization properties of the 
poly(dG)-poly(dC), poly(dA)-poly(dT) DNA polymers 
and the mixed cases. Figure 5(a) and (b) show the local- 
ization length and the Lyapunov exponent in the three 
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FIG. 4: (a)The on-site energy EncV and (b)transfer integral 
Vnn+ieV as a function of the base pair site n. The parametric 
plot E„ versus Vnn+i is shown in (c). W — 0.1 and the other 
parameters are given in text. The unit of the energy and 
the spatial length are eV and the number of nucleotide base 
pair(6p), respectively, throughout the present paper. 
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FIG. 5: Comparison (a)localization length and (b)Lyapunov 
exponent in poly(dG)-poly(dC), poly(dA)-poly(dT) and the 
mixed DNA polymers. = 0.1, So = 36° and iV = 2^^ 
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FIG. 6: The localization length as a function of the energy 
for several system sizes iV(= 2^*^, 2", 2^^) in the (a)poly(dA)- 
poly(dT) DNA polymers, (b)poly(dG)-poly(dC) DNA poly- 
mers. W = 0.1 and 9o = 36°. 



tic conductance of DNA in water solutions, which is also 
temperature-independent. Zalinge et al. have tried to 
explain the latter effect, using a kind of acoustic phonon 
motions in DNA duplexes, which seems to be plausible, 
but not the only possible physical reason [fl |45|. We 
will propose an alternative explanation for the observed 
temperature-independent conductance, based upon our 
numerical results. 

Generally, in the case of quantum diffusion the tempo- 
ral evolution of the electronic state vector | $) is described 
by the time-dependent Schrodinger equation ih^^^ = 
Hei{t)\^), which then becomes 



Vnn + l{t)4>'n.+ l - K-ln (i)0)(J-7) 



types of polymers with W = OA. In low energy regions, 
the localization length in the poly(dG)-poly(dC) DNA 
polymer is larger than that in the poly(dA)-poly(dT) 
DNA polymer. The system size dependence of the lo- 
calization length is given in Fig. 6 in relation to the reso- 
nance energy. Indeed, the localization length of the DNA 
polymers is larger than ^ > 2000 [fop] in almost all the 
energy bands for all models; it is much larger than the 
system size of the oligomer used in the experiments. As 
is seen in Fig. 6, the smaller the size of the system the 
more complex the resonance peaks become. 



D. Quantum Diffusion in the Fluctuating 
Environment 



In this subsection we numerically investigate quan- 
tum diffusion of electrons in the Hennig model of 
poly(dG)-poly(dC) and poly(dA)-poly(dT) with a dy- 
namical disorder [l^. We assume that the diffusion is 
caused by a colored noise associated with the stochas- 
tic dynamics of distances r(t) between two Watson-Crick 
base pair partners: (r(t)r(t')) = roexp(— |i — t'l/r). 
These fluctuations can be regarded as a stochastic pro- 
cess at high temperature, with phonon modes being ran- 
domly excited [42l. l43j. In the model, the characteristic 
decay time t of correlation can control the spread of 
the electronic wavepacket. Interestingly, the white-noise 
limit r — >■ can in effect correspond to a sort of mo- 
tional narrowing regime(see, examples), because we find 
that such a regime causes ballistic propagation of the 
wavepacket through homogeneous DNA duplexes. Still, 
in the adiabatic limit r — > oo, DNA electronic states 
should be strongly affected by localization. The ampli- 
tude ro of random fiuctuations within the base pairs (the 
fiuctuation in the distance between two bases in a base 
pair) and the correlation time r, are very critical param- 
eters for the diffusive properties of the wavepackets. It 
is interesting to find the ballistic behavior in the white 
noise limit vs the localization in the adiabatic limit, since 
there is a number of experimental works observing ballis- 



where 4>ri = and the effective Planck constant 

heff = 0.53. We redefined the scaled dimensionless vari- 
ables En{t) and Vnn+i{t) in Eq.(l) such that -p^ — > -E„, 

— Vnn+i. We used mainly the 4th order Runge- 
Kutta-Gill method in the numerical simulation for the 
time evolution with time step = 0.01. In some cases, 
we confirmed ourselves that the accuracy in the obtained 
results is in accord with the one in the results that are 
gained by the help of the 6th order symplectic integrator 
that is the higher order unitary integration. 
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FIG. 7: Diffusion rate D{t) as a function of for several 
fluctuation strengths rno — 1, 3, 5 at (a)A-T model and (b)G- 
C model, respectively 
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FIG. 8: Short-time behavior of m{t) of A-T, G-C, and mixed 
models with r„o = 1-0 at (a)r — 1 and (b)r — 0.01. 



It has been demonstrated that the motional narrow- 
ing affects the localization in the poly(dG)-poly(dC) and 
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poly(dA)-poly(dT) DNA polymers. In either case of the 
model DNA polymers, the temperature-dependence be- 
comes virtually suppressed when the motion of the wave 
packet is characterized by the ballistic propagation. 

We have also investigated the temporal diffusion rate 
in almost all the diffusive ranges. (See Fig. 7.) We found 
that the diffusion rate of the A-T model is larger than 
that of the G-C model at comparatively low tempera- 
tures. Interestingly enough, for relatively high tempera- 
tures in the diffusive range of the wavepacket motion the 
difference between the two DNA systems gets smaller. 

Figure 8 shows the short-time behavior in the cases 
of T = 1 and T = 0.01 depicted on a larger time scale. 
It follows that in the short-time behavior the (dG)i5 — 
(dC)i5 case is more diffusive than its {dA)i5 — {dT)i5 
counterpart within the range from which the spread of 
the wavepacket is y/m ^ 15. 

We used periodic sequences, which means that Eq 
and Vq are constant values as the static parts of the 
on-site and hopping terms for poly(dG)-poly(dC) and 
poly(dA)-poly(dT) DNA polymers, respectively. This 
includes the mixed model as well. Then the motional 
narrowing for dynamical disorder makes time-evolution 
of the wavepacket ballistic. However, we remark that the 
motional narrowing strongly localizes the wavepacket if 
we use a disordered sequence for the static parts of £"„ 
and/or Vnn+i- 



simple model to investigate electronic diffusion coupled 
with lattice vibrations. 

In Fig.9(a), the MSD is shown for the poly(dG)- 
poly(dC), the poly(dA)-poly(dT) and the mixed DNA 
polymers. We find that all the cases exhibit the normal 
diffusion of electron without any stochastic perturbation. 
As shown in Fig.9(b), the diffusion rate decreases as the 
perturbation strength e increases. As a result the coupled 
motion of charges and the lattice breathers connected 
with the localized structural vibrations may contribute 
to the highly efficient long-range conductivity. (See ap- 
pendix D.) 




FIG. 9: (a) MSD m(t) in A-T, G-C, and mixed models per- 
turbed by M = l,e = 0.5. (b) MSD m{t) of A-T model 
perturbed by M = 1 for e = 0.5, 0.7, 1.0, 2.0 from above. We 
set r„o = 1. 



E. Quantum Diffusion Coupled with Vibrational 
Modes 

Let us denote r„ (t) by the stretching vibration of W- 
C bonds at the n site as before. The external harmonic 
perturbation is equivalent to the coupling with quantum 
linear oscillators or with phonon modes in solid state 
physics m, lU . We replace the disordered fluctuation 
r„(t) with the harmonic time-dependent one such as 

M 

rn{t) = r„o + ^"^)' (18) 

1=1 

where M and ti are the number of the frequency compo- 
nent and the strength of the perturbation, respectively. 
The initial phases {9n^} e [0, 27r] at each site n are ran- 
domly chosen. In the numerical calculation of this sec- 
tion, we take — -^^j, for simplicity, and take incom- 
mensurate numbers of ^ 0{1) as the frequency set. 
In the limit of M — > oo, the time-dependent perturba- 
tion approaches the stochastic fluctuation as discussed in 
the last subsection. In particular, we can regard the ap- 
proximation as an electronic system coupled with highly 
excited quantum harmonic oscillators. One of the ad- 
vantage of this model is that although the number of au- 
tonomous modes is limited due to the computer power, 
we can freely control the number M of frequency com- 
ponents in the harmonic perturbation. This is also a 



V. TIGHT-BINDING MODELS FOR THE 
LADDER SYSTEMS 

In this section, to investigate the energy band struc- 
ture for periodic sequences and the localization proper- 
ties for the correlated disordered sequences, we intro- 
duce the ladder models of DNA polymers. Although in 
Refs.fTTI, [Tsl . [2l| . we assumed the ladder structure con- 
sisting of sugar and phosphate groups, the model can 
be applicable to the DNA-like substances such as the 
ladder structure of bases pairs without sugar-phosphate 
backbones. When we apply the model for the sugar and 
phosphate chains, the chain A and B are constructed 
by the repetition of the sugar-phosphate sites, and the 
inter-chain hopping Vn at the sugar sites come from the 
nucleotide base-pairs, i.e., ^ — T or G — C pairs. (See 
Fig.lO(a).) 



A. The Ladder Models and the Dispersions 

Since there is only one tt— orbital per site, there are 
totally 2N tt— orbitals in the ladder model of the DNA 
double chain. Let us denote by (pnifpn ) the tt— orbital 
at site n in the chain A{B). By superposition of the 
TT— orbitals, the Schrodinger equation H\'^) = E\'^) be- 
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FIG. 10: Models of the double strand of DNA. (a) the two- 
chain model, (b) the three-chain model we adapted in the 
main text. 



comes, 

Bn+l,n4'n+l + Bn.n-l4>n-l + Bn,n't'n + K0„ = E4>^ , 



where An+i^n {B„+i^n) means the hopping integral be- 
tween the nth and (n + l)th sites and A„ „ {Bn,n) the 
on-site energy at site n in chain A (B), and Vn is the hop- 
ping integral from chain A{B) to chain B{A) at site n, re- 
spectively. Furthermore it can be rewritten in the matrix 
form, $„+i = MA, where = (0;^, <_i)t. 
The transfer explicit matrix is given in appendix B. Gen- 
erally speaking, we would like to investigate the asymp- 
totic behavior {N — ^ oo) of the products of the matrices 
Md=2 (iV) = ^LiMk = MnMn-i • • • Ml. 

We consider the band structure for the periodic case by 
setting An+i^n = Bn+i^n = ci{b) at odd (even) site n and 
Vn = v{0) at odd(even) sites, Ann = Bnn = a{Ann = 
Bnn = P) for odd(even) site n, for simplicity. 

A simple way to solve the above equations is to use the 
Bloch theorem: 



A2ks 



(19) 



where s is the length between the adjacent base groups 
and is assumed to be equivalent to the length between 
the adjacent sugar-phosphate groups and the wavenum- 
ber k is defined as — t?^ < k < (We take s = 1 

Zs — — 2s ^ 

in the following calculations.) The Schrodinger equation 
becomes 



M 



M^n, 
( ^ 



a + be-''') 





-{a + be''') 
E -a 


—V 






-{a + be-'''] 





~v 

-{a + be''') 
E -a 



(20) 
(21) 



where $„ = (04+i> </'2l; '/'fn)*- Denote the deter- 

minant of M by D{M). Solving the equation D{M) = 
for E', we obtain the energy dispersion of the system. 
Here we show the simple case of a = /3 = 0, 

E^^\k) = i [v ± ^ + 4(a2 + 62) + 8a6 cos 2fc) ,(22) 

E^^\k) = i (-i>± V-y^ +4(a2 + 62) +8a6cos2fc) (23) 

E^\k) ( E'^\k) ) stand for the upper(lower) bands 
in channel ±, respectively [ll|. The lowest and upper 
middle (the lower middle and highest) bands correspond 
to bonding (antibonding) states between the adjoint or- 
bitals in the interchains. The more detail of the energy 
band for general case is given in appendix C. In Fig. 11 (a) 
the energy band structure is given with varying the in- 
terchain hopping V. Figure 11(b) shows the cross-section 
view at V = 1. There is a band gap Eg(y) at the cen- 
ter in between the lower and upper middle bands in the 
spectrum for the whole range of v when a — (3 — Q. 

Eg{v) = E^+\n/2) - E^~\7t/2) = V4(a - 6)2 + (24) 



The other band gap Ag{v), 



Ag{v) = v + - (v/4(a - 6)2 + v^- ^A{a + bf + v^) 



appears in between the lowest and the lower middle bands 
(the upper middle and highest bands) when v > Vc = 
2ab/^Ja? + b^ (otherwise, it is negative and therefore 
semimetallic). There is a transition from semiconduc- 
tor to semimetal as the tt— electron hopping between the 
nitrogenous bases of nucleotide is increased. 

The two-chain model can be easily extended to the 
three-chain case. (See Fig. 10(b).) C„+i,„ means the hop- 
ping integral between the nth and {n + l)th sites and 
Cn.n the on-site energy at site n in chain C, and Vn and 
Un are the hopping integral between the chains. The 
Schrodinger equation H\'^) — E\'^) becomes. 
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f Un<t>n 
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It can be also rewritten in the matrix form, $„+i = 
Md=3(n)$„ where = (/)^, (/>^_)* 

and the expHcit transfer matrix is given in the appendix 
C. 

For simphcity, we set An+i,n = Bn+i,n — ol and 
Cn+i,n = c at site n, and Vn = Un- We can analyti- 
cally obtain the energy band structure for the three chain 
model by using the Bloch theorem. 

= e^'^n . = e'"=^</.^>,f+i = e^'^^^/.^f . (27) 



B. Localization in the Ladder Models 

As was discussed in Sect. 2 and Sect. 3, the sequences 
of the realistic DNA polymers are not periodic and ac- 
company a variety of disorder. Generally, the random- 
ness makes the electronic states localized because of the 
quasi-one-dimensional system and affects electronic con- 
duction and optical properties, and so on. In the present 
section, we investigate the correlation effect on the local- 
ization properties of the one electronic states in the dis- 
ordered, two-chain (ladder) and three-chain models with 
a long-range structural correlation as a simple model for 
the electronic property in the DNA[i3|. The relationship 
between the correlation length in the DNA sequences and 
the evolutionary process is suggested. Moreover, it is 
interesting if the localization property is related to the 
evolutionary process. 

Figure 12 shows the DOS as a function of energy for 
the binary disordered systems. The sequence of the in- 
terchain hopping Vn takes an alternative value Wgc or 
Wat- We find that some gaps observed in the periodic 
cases close due to the disorder corresponding to the base- 
pair sequences. 

Figure 13(a) shows the energy dependence of the Lya- 
punov exponents (71 and 72) for some cases in the asym- 
metric modified Bernoulli system characterized by the 
bifurcation parameters Bq and Bi controlling the corre- 
lation. (See appendix G for modified Bernoulli system.) 
They are named as follows: Case(i):i3o = 1.0, i?i = 1.0, 
Case(ii):Bo = 1.0, Bi = 1.9 and Case(iii):Bo = 1-7, Bi = 
1.9. Apparently, the case(i) is more localized than the 
cases (ii) and (iii) in the vicinity of the band center 
\E\ < 1. The comparison between the case(ii) and the 
case (iii) shows the effect of asymmetry of the map on the 
localization. The ratio of the G-C pairs Roc ^ 0.2 for 
the case(ii), Rgc 0.47 for the case(iii). In the energy 
regime \E\ > 1, the Lyapunov exponent 72 in the case(ii) 
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FIG. 11: Energy bands of the decorated ladder model, (a) 
The energy bands as a function of V. (b) The snap shot of 
the energy bands when V — 1. k means the wave vector in 
units of - such that —1.0 < k < 1.0, V means the 7r-electron 
hopping integral between the inter chain sites, and E means 
the energy in units of = 1. Here we have taken the values 
a = 0, /3 = 0, a = 0.9, 6 =1.2. 



is smaller than the one in the case(iii) in spite of the same 
correlation strength Bi. As a result, we find that in the 
DNA ladder model, correlation and asymmetry enhance 
the localization length ^(= I/72 ) of the electronic states 
around \E\ < 1, although the largest Lyapunov expo- 
nents 7i do not change effectively at all. Figs. 13(b), 
(c) and (d) show the nonnegative Lyapunov exponents 
in the real DNA sequences of (b) IICh-22, (c) bacterio- 
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phages of E. Coli and (d) histon proteins. In the case of 
HCh-22, we used two sequences with N = 105, extracted 
from the original large DNA sequences. The result shows 
that the Lyapunov exponents do not depend the detail of 
the difference in the sequences of HCh-22. Although the 
weak long-range correlation has been observed in HCh- 
22 as mentioned in the introduction, it does not affect on 
the property of localization. The sequences we used are 
almost symmetric (Rgc ~ 0.5). 
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FIG. 12: DOS as a function of energy for the binary disor- 
dered cases in the two-chain model. (a)WAT = 2.0, Woe = 
1.0, a = 1.0, b = 1.0. (b) Wat = 2.0, Wgc = 1.0, a = 1.0, 6 = 
0.5. The on-site energy is set st Ann = S„„ — 0. 



In addition, in the numerical calculation we set the on- 
site energy as Ann = Bnn = Gnn = for simplicity. The 
sequence {C„„+i} can be also generated by correspond- 
ing to the base-pairs sequence. The localization prop- 
erties in the simple few-chain models with the on-site 
disorder have been extensively investigated (481] ■ Figure 
14 shows the energy dependence of the Lyapunov expo- 
nents in the three chain cases. We can observe that all the 
Lyapunov exponents 7^(1 < d) are changed by the corre- 
lation. The least nonnegative Lyapunov exponent 73 is 
diminished by the correlation. In particular, it is found 
that the correlation of the sequence enhances the localiza- 
tion length defined by using the least non-negative Lya- 
punov exponent. We can see that the localization length 
diverges at the band center E = Q. 

In Fig. 15(a) the result for asymmetric modified 
Bernoulli system is shown. Apparently, the correlation 
and/or asymmetry of the sequence effect a change in the 
second and third Lyapunov exponent. In contrast, al- 
though the global feature of 71 is almost unchanged, the 
local structure of the energy dependence is changed by 
the change of Bq. Fig. 15(b) shows the results in phage-fd 
and phage-186 in the three-chain model. It is found that 
the structure of the energy dependence around < 2 
is different from that in the artificial sequence by the 
modified Bernoulli map. 

The charge transfer efficiency based on mismatch and 
correlation effect for some genomic and synthetic se- 
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FIG. 13: Lyapunov exponents (71,72) as a function of en- 
ergy in the ladder model, (a) Modified Bernoulli model, where 
cases(l):Bo = 1.0, Bi = 1.0, Case(2):Bo = 1.0, Bi = 1.9 and 
Case(3):Bo = 1.7, Bi = 1.9 are shown. (b)HC-22. (c) bac- 
teriophages of E. coli (phage-A, phage-186). (d) early his- 
tone HI and late histone HI. The on-site energy is set at 
Ann = Bnn = 0, fl = —1.0, h = —0.5. The size of the se- 
quence is TV = 10^ for (a), = 10^ for (b), iV = 48510 
for the phage-A in (c), N — 30624 for the phage-186 in (c), 
iV = 787 for the early histone HI in (d), and N = 1182 for 
the late histone HI in (d). 




FIG. 14: Lyapunov exponents 7^(1 = 1,2,3) as a function of 
energy in the correlated three-chain model. The parameters 
Wat = 2.0, Wat-at = 1.0, a = 1.0,6 = 0.5. The on-site 
energy is set at Ann = Bnn = C„„ = 0. 



quences are also investigated [4| 
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FIG. 15: Lyapunov exponents (71, 72, 73) as a function of en- 
ergy in the three-chain model, (a) Modified Bernoulli model, 
(b) bacteriophages of E. coli (phage-/d, phage-186). The pa- 
rameters are same as ones in Fig. 14 except for on-site energies 
of C-chain (Cn,n — 0). 



C. Polaron Models for a Double Strand of DNA 

The overlap of the electronic orbitals along the stacked 
base pairs provides a pathway for charge propagation 
over 50A and the reaction rate between electron donors 



and acceptors does not decay exponentially with dis- 
tance. A multiple-step hopping mechanism was recently 
proposed to explain the long-range charge transfer be- 
havior in DNA. In this theory, the single G-C base pair 
is considered as a hole donor due to its low ionization 
potential if compared to the one on the A-T base pairs. 
Besides the multiple-step hopping mechanism, polaron 
motion has also been considered as a possible mecha- 
nism to explain the phenomena of DNA charge transfer. 
The charge coupling with the DNA structural deforma- 
tions may create a polaron and cause a localized state. 
The polaron behaves as a Brownian particle such that 
it collides with low energy excitations of its environment 
that acts as a heat bath. For the most physical systems, 
the acoustical and optical phonons are the main lattice 
excitations as in the single-chain case given in Appendix 
D. 

In this section, we present the formation of localized 
polarons due to the coupling between charge carriers and 
lattice vibrations of a double strand of DNA(i^. These 
polarons in DNA act as donors and acceptors, which ex- 
hibit an extrinsic semiconductor character of DNA. The 
results are discussed in the context of experimental ob- 
servations. 

Let Xn and ?/„ be configuration angles for rotation of 
the nth base group along the course of the backbone 
chains A and B of DNA ladder, respectively. For the 
sake of simplicity, we assume that all masses Mn of the 
nucleotide groups and the bond lengths dn between S and 
Bn groups are identical so that Af„ = M and d„ — d, giv- 
ing the moment of inertia /„ = Iq = M d^ . In this case, 
the lattice Hamiltonian Hph can be given by 



A. 



K , 



(28) 



where A, K and S are the parameters for the rotational 
energy, the stacking energy and the bonding energy of 
the bases, respectively. We note here that various gen- 
eralizations and modifications of the model are straight- 
forwardly possible by adapting a different combination of 
the interactions and the hopping integrals. In the follow- 



ing, we give result only for HOMO. We can obtain the 
result for LUMO in simple replacements. 

On the other hand, generalizing the idea of Holstein for 
the one-dimensional molecular crystal to our case of the 
double strand of DNA, the tight bingding Hamiltonian 
for electrons is given by 
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N 



n=l 
N 



For the case of LUMO, we can obtain the resuh hy a b 
in Here a^^b^^) is the electron creation operator 
in the HOMO (LUMO) at the nth nucleotide group of 
P, S and Bn for the chain A. e^(a;„) and e^(j/n) are the 
electron-lattice coupling Hei-ph given by 



4{y) = 4- FhV, 



(29) 
(30) 



respectively. We can also obtain e;^(x„), e£(j/„) for 
LUMO. Moreover, we assume the relation: 



(31) 



since this condition is realized in most of the DNA sys- 
tems and guarantees the semiconductivity of the DNA. 
The total Hamiltonian Htot is given by 



Htot = Hph + H^i + H^i, (32) 
while the wavefunction |$) given by 

N 

s=A,B n=l 

where al^„ (6jj) means the electron creation operator at 
site n in chain A (B). 4>n{4'n) represents the HOMO 
and fni'Pn) the LUMO, at the nth nucleotide site in the 
chain A(B), respectively. Applying to the Schrodinger 
equation Htot\^) = Etot\^) we obtain the following equa- 
tions for the HOMO states of nucleotide groups in the 
double strand of DNA: 



D. DC-Conductivity of the Double Strand of DNA 

We exclusively focused on the low-energy transport, 
when the charge injection energies are small compared 
with the molecular bandgap of the isolated molecule 
which is of the order of 2 — 3eV. In the experiment of the 
conductance property of the DNA, temperature depen- 
dence is important. Finite temperature can also reduce 
the effective system size and leads to the changes in the 
transport property. Moreover, the effects of the stacking 
energy and of temperature can be considered by intro- 
ducing the fluctuation in the hopping energy such as the 
Su-Schriefer-Heegar model for polyacetilene [50l[5l|. 

We consider DC-conductivity of periodic DNA se- 
quence based on the LUMO, HOMO band of tt— electrons 
and small polaron [l^. It is known that DNA behaves 
n(p)— type extrinsic semiconductor, where the donors 
(acceptors) for the LUMO (HOMO) band are positively 
(negatively) charged small polarons with the total num- 
ber Nd{Na) and the energy ea = E'/ ^ - FI/IqujI 

— E^ = + Ffi/IoujQ, where E^{Ec) is the valence 
(conduction) band edge. Following the standard argu- 
ment, denote by ric and n/^ the numbers of electrons in 
?iie conduction (LUMO) band and in the donor levels, 
respectively; and denote by py and pA the numbers of 
electrons in the valence (HOMO) band and in the ac- 
ceptor levels, respectively. We now have the relation: 

TT-c + Ud = Nd - Na + Pv + Pa, whcrC 



N, 



rid 



(«d-M) -f 1 



5/3(M-«d) + 1 ' 



(36) 



n+l 



^l} + 



FHXn<, 



where f3 = 1/KbT, fj, is chemical potential of the system. 
If we suppose ed — fJ^ >> ksT, ^ — €a >> ksT and 
Ud « Nd, Pa « Na, then 



where E = Etot~ Hph{xn,Vn), Vn = v{\x„-y„\). For the 
LUMO states of nucleotide groups in the double strand 
of DNA, the similar equations are given by replacing 
and H by ip and L, respectively. Moreover, according to 
argument of Holstcin, we find the following coupled dis- 
cretized nonlinear Schrodinger equation(DNSE) describ- 
ing the electronic behavior under the lattice vibrations 
in the ladder model of DNA. In Appendix E, we give the 
derivation of the coupled DNSEs. The single-chain ver- 
sion of the DNSE and some comments on the physical 
meanings are also given in Appendix D. 



where Ui = y/N^erl^^f/^, fn 
\og{Py/N,)/2(3. 

The DC conductivity is given 



e^ncTe 
a = (= 



e PvTh 
rrih 



(37) 

{E, + Ey)l2 + 



(38) 



by Drude formula. This suggests that if Tg — const then 
the temperature dependence of a comes from ric, while 
if n^ = const then it comes from t^. Since the temper- 
ature dependence of Tg comes from other mechanisms of 
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scattering such as the activation of polaron motion con- 
sidered by Yoo et al., we can assume = Toe~^^'^ where 
Ea is the activation energy. Hence, we have 

a = exp —phja (39) 

me 

If we adopt the logrr vs 1/T plots, wc found the strong 
temperature dependence found by Tran et al. and Yoo 
et al. (See Fig.l6.) 




FIG. 16: The DC conductivity is shown for the n-type (or 

p-type) extrinsic semiconductor where the energy gap Eg = 
0.33eF, the hopping energy for the HOMO (or LUMO) band 
t = 0.2eV and the activation energy for the polaron hopping 
of (a) Ea = O.OOli, (b) Ea = O.OSi and (c) Ea = 0.05t are 
taken, respectively. Nc = Pv = 10* and n = No - Na = 10^ 
are assumed. 
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FIG. 17: Energy gap Eg{T) of the n— type (or p—type) ex- 
trinsic semiconductor is shown for electron conduction (or 
hole conduction) where the energy gap Eg = 2eV, t = O.leV, 
Nc = Pv = 10® and n = 10® are used, respectively. 

It is found that the band gap is reduced by the for- 
mation of the double strand of DNA and small polarons 



exist to behave as localized donors and acceptors in the 
DNA double helix. (See Fig. 17.) This extrinsic semi- 
conducting nature of DNA qualitatively explains many 
experimental results that are observed recently. 

As was stated in the introduction, it is very well-known 
that in solid state physics, the extrinsic semiconductor 
character is added by impurity atoms that are exerted 
from outside into host materials. For example, if the 
host material is a homogeneous crystal of Si or Ge with 
four valence bonds, then the impurities are As atoms 
with five valence bonds and boron(B) atoms with three 
valence bonds, which then provides an inhomogeneous 
material with impurity levels. The As atoms play a role 
of donors with positive charge such that the system be- 
comes an n-type semiconductor, while the B atoms play 
a role of acceptors with negative charge such that the 
system becoms a p-type semiconductor. Such the system 
of n- or p-type semiconductor exhibits an extrinsic semi- 
conductor character. On the other hand, in our systems 
of DNA semiconductors, there exist no such impurities 
exerted from outside; but already there exist inside the 
DNA a complicated arrangement of bases of A, G, C and 
T, such that among the bases, each one of bases may re- 
gard other bases as impurities. Hence, a kind of extrinsic 
semiconductive nature may appear as a result of self or- 
ganization of the base arrangement. This is the meanig 
of our terminology of "extrinsic semiconductor" for the 
DNA systems. Therefore, it is more preferable for us to 
put " self-organized" in front of " extrinsic semiconductor" 
for DNA such as "self-organized extrinsic semiconduc- 
tor" in order to represent the semiconducting character 
of DNA. 

We would like to note that in the rf— dimensional dis- 
ordered systems, the activated hopping between local- 
ized states, i.e., the variable range hopping, could be a 
dominant mechanism for the dc— conductivity, and the 
temperature dependence is governed by 

a = c7oe-(^)"^, (40) 

where Tq = SWaj/kB- Here W is the energy difference 
between the two sites, 7 is Lyapunov exponent, and a the 
distance between the nearest neighboring bases. It is rea- 
sonable to expect that the experiments of charge trans- 
port in A— DNA suggest such temperature-dependence 
a oc e~^f^^^ at relatively low temperature due to the 
Id aperiodic base-pair sequence. The main contribu- 
tion was given by the interaction with water molecules 
and not with counterions. Further, polaron formation 
was not hindered by the charge- solvent coupling. And 
the interaction rather increased the binding energy (self- 
localization) of the polaron by around 0.5eV, which is 
much larger than relevant temperature scales. 

Before closing this section, we give brief comments con- 
cerning the measurements of electronic conductivity in 
DNA again. The DNA is quasi-lD systems coupled with 
the environment including the surrounding substrate and 
contact leads. In general, Landauer-Buttiker formula us- 
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ing transfer matrix method and Kubo-Geenwood formula 
using linear-response theory have been used to estimate 
the transport properties in the quantum systems. The 
transfer matrix method is straightforward for the quasi- 
ID systems, however, it is inconvenient for the compli- 
cated cases involving the environmental effects although 
the two formulae are equivalent at least for single parti- 
cle cases 52, 5^]. There is an indication that the transfer 
matrix method is also effective even for law transmission 
coefficient due to the contact effect, despite a very good 
charge transfer along the DNA sequence j53 |. 



VI. CONDUCTION AND PROTON TRANSFER 

Gene mutations sometimes cause human disorders such 
as cancer. Simultaneously, gene mutation can be a driv- 
ing force for processes of biodiversity and evolution j55|. 
When a system contains hydrogen bonds, proton mo- 
tion also needs to be considered. Instead of oscillatory 
motions in a single-well potential, protons can tunnel 
from one side to another in a double-well potential of 
the hydrogen bond. This proton tunneling causes an 
interstrand charge hopping. And it could generate the 
tautomeric base pairs and destroy the fidelity of the W- 
C base pairs. 

Lowdin proposed that proton tunneling contributes to 
the formation of rare tautomeric of DNA W-C base pairs 
whose accumulation could result in DNA damages, point 
mutations, and even tumor growth [56l|. The argument re- 
lies on the assumption that the rare tautomers are more 
stable, and once the intermolecular proton transfer oc- 
curs, lifetime of the tautomers is significantly larger than 
time of DNA repair. Indeed, the results from the quan- 
tum chemical and statistical mechanical calculations in- 
dicate that at room temperature at least the GC tau- 
tomers (G* = C*) would have a sufficient lifetime to 
cause the DNA damage through mismatches of the W-C 
base pairs. 

Ladik speculated that semiconductivity of DNA might 
be related to the origin of cancer due to transmutation 
of genes 0. Very recently, Shi et al. have also reported 
correlation between the cancereous mutation hot spots 
and the charge transport along DNA sequences [s^l- 

The proton transfer causes fluctuation of the poten- 
tial and more or less affects on the localization/transport 
property of the charged carriers. At the same time, the 
proton transfer suffers from the lattice vibration and the 
electron transport along the base sequences of DNA. Ac- 
cordingly, we should treat the coupling between the pro- 
ton transfer and the phononic and electronic states when 
we investigate the stability of the tutomatic states(the 
excited states). Recent theoretical studies have shown 
that charged protonated base pairs display smaller acti- 
vation barriers, which make the proton transfer and the 
tunneling more facile. 

Chang et alJ^] have given the effective ID model 
Hamiltonian Htot = Hei-ph(t) + + Hei-a, where 



Ha^ = J2n Han represents Hamiltonian describing the 
proton transfers. Here Hamiltonian Hei-ph includes the 
electron and phonon modes and the coupling. For the 
proton motion in hydrogen bonds, they used a two-level 
system in order to describe the tunneling behavior at 
molecular site n 

Han^\{~€pK+tpan), (41) 

where cr^ and cr^ are Pauli matrices assigned at site n, Cp 
the energy bias between the two localized proton states 
and tp the tunneling matrix element between the states. 
They have modeled the coupling between the protons in 
the hydrogen bond and the charges in the DNA strand 
as 

H,l_a^Y.^an{'yl-l)ClCn (42) 
n 

where is the coupling intensity. When the proton is 
in the lower energy state or there is no charge around the 
hydrogen bond ((Cj,C„) = 0), the coupling vanishes. 

Indeed, the normal G-C base pair is in the lower energy 
state and its tautomeric form is in an excited state. When 
Cp » tp, the probability of having a tautomeric form 
is extremely small. However, the cation of a G-C base 
pair has almost the same energy as that of its tautomeric 
form ep < tp. In this case, the proton state becomes 
delocalized in hydrogen bonds. 

The symmetry of potential well and the hight of barrier 
essentially affects on the tunneling probability between 
the potential wells. The probabilistic amplitude of pro- 
ton could influence the transport of the charged carriers 
though the Coulomb interaction. In addition, the two- 
level approximation will be broken down if the chaotic 
motion occurs in the dynamics due to the coupling with 
lattice vibrations (ssj . 



VII. SUMMARY AND DISCUSSION 

We briefly reviewed our recent works with concern- 
ing the electronic states and the conduction/transfer in 
DNA polymers. In Sect. 3, based on the Hitckel model, 
we have discussed the basics of quantum chemistry where 
the electronic states of atoms in biochemistry such as 
C, N, O, P and the electronic configurations of the ni- 
trogenous bases, the sugar-phosphate groups, the nu- 
cleotides and the nucleotide bases are summarized, re- 
spectively. In Sect. 4, based on the tight-binding approx- 
imation of the DNA sequences, we have investigated the 
localization properties of electronic states and quantum 
diffusion, using a stochastic-bond-vibration approach for 
the poly(dG)-poly(dC), the poly(dA)-poly(dT) and the 
mixed DNA polymers within the frameworks of the po- 
laron models. At that time, we assumed that the dynam- 
ical disorder is caused by the DNA vibrational modes, 
which are caused by a noise associated with stochastic 
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dynamics of the distances r{t) between the two Watson- 
Crick base pair partners. In Sect. 5, we have mainly in- 
vestigated the energy band structure and the locahza- 
tion property of electrons in the disordered two- and 
three-chain systems (ladder model) with long-range cor- 
relation as a simple model for electronic property in a 
double strand of DNA by using the tight-binding model. 
In addition, we investigated the localization properties of 
electronic states in several actual DNA sequences such as 
bacteriophages of Escherichia coli, human chromosome 
22, compared with those of the artificial disordered se- 
quences. In Sect. 2 and 4, we gave a brief explanation 
for DNA bases sequences and gene mutations, respec- 
tively, which are related to the DNA carrier conduction 
phenomena [59-'62']. The relationship between the long- 
range correlations and the coherent charge transfer in the 
substitutional DNA sequences has been studied, based on 
the transfer matrix approaches [63[ [67| . 

Recently, the tight-binding models for the DNA poly- 
mers have been extended to the decorated ladder models 
and the damaged DNA models by some groups. Further- 
more, the I-V characteristics of the ladder models cou- 
pled with environments has been investigated as simple 
models for DNA polymers |63 - I66j . 



Appendix A: Hiickel Parameters 

Sandorfy's form ula [28|] . Streitwieser's formula [30|. and 
Mulliken's formula[2^ are given by following relations. 



XX - Xc 
Xc 



X 4.1/3, 



"x 



5, 



cx 



f-cx 



Sec 



(Al) 
(A2) 
(A3) 



respectively. See Ref.[5| for the details. 



Appendix B: Transfer Matrix Method 

Let us define the four-dimensional column vector in 
the main text as <&„ = ((/)J^, t/)^, (/)^„j)*. Then, it 

can be rewritten in the following form: 
where 



M„ ]■ (Bl) 



M„ is the 4x4 transfer matrix with the 2x2 matrices: 



According to the sequence of A'' segments, we have to take 
the matrix product M{N) = MnMn-i ■ ■ ■ Afi, which is 
also a 4 X 4 matrix. When the double chain system is 
periodic, we adopt the Bloch theorem: 

0n+iV = P'/'n , €+N = (B4) 

where p = e^^-^ and k is the wavevector. Then, we find a 
4x4 determinant D{p) that is a fourth order polynomial 
of p: 

D{p) =det[M{N)^ ph] 
= p^ ~ aip^ + a2p^ — a^p -|- 04 = 0. (B5) 

It provides the wave vector k in the system such that 
k = 2TTj/N for j = -N/2, ■ ■ ■ ,N/2. Here, if four roots 
are written as pi, p2, ps, and p4, then the following con- 



ditions hold 

4 4 
ai = trM{N) = ^ p^, 02 = ^ PiPj, 

i—l i<j — l 

4 

0-3= ^ piPjPk,a4, = det M{N) = piP2P3Pi- (B6) 

i<j<fe=l 

Using a physical intuition, if an electron propagation 
with k along one direction in the double chain is rep- 
resented by p, then the reverse propagation with — fc is 
represented by Therefore, the latter should be also 
accessible, since the choice of the direction of the coordi- 
nate system is arbitrary. Hence, p~^ must be an eigen- 
value of D{p) =0 such that 

D{p-^) = p-^ia^p"- ~ asp^ + a2P^ - aip + I) = 0.(B7) 
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This situation imposes the particular condition on the 
matrix M: 



mUm = J, 



-1 

1 



(B8) 



(B9) 



where M'^ means the Hermitian conjugate of M and 
is the 2x2 zero matrix, respectively. This property is 
called the symplectic structure of M, where we have 



D{p) = p^D{p-'), 



(BIO) 



from which we find ai = as, 04 = 1. Thus, M belongs to 
5L(4,R). By using this property and dividing D{p) by 
p^, the biquadratic equation is reduced to the quadratic 
equation: 



X — aix + a2 — 2 = 0,x = p + 



1 



(Bll) 



where £>_ (E) [D+ (E)] means the DOS contributed from 
the bonding (antibonding) channel — (+), respectively. 
It agrees with the result on the tight-binding model for 
the ladder structure. Physically speaking, the — (+) 
channel means the bonding (antibonding) states between 
two parallel strands of the DNA. 

Appendix C: Transfer Matrices for Ladder Systems 

In this appendix, we give the explicit expression of 
the energy band for the decolated ladder model given in 
Fig. 10(a). In the unit cell, the period is taken as = 2 
and it contains four tt— orbitals. We apply the result in 
Appendix B for the case. Let the transfer matrix method 
be 



M 



A V 
U B 



(CI) 



M is a 4 X 4 transfer matrix with 2x2 matrices: 



Therefore, its two roots are given as 



x± 



\ (ai ± Vd) , D^aj-4a2 + 8. (B12) 



Since from Eq. (jB6P and tr{AP) — J2i=iPh ^'^ obtain 
aj - 2a2 = tr{M'^). 

Now we can state a simple scheme to obtain the spec- 
trum: If an energy E satisfies 



x± = 2 cos kN, 



(B13) 



then the energy is allowed, otherwise it is forbidden in 
channel ±, respectively. This is a generalized version of 
the Bloch condition for the single linear chain system 
with the 2x2 transfer matrix M where 



{E~P)iE^a) _ a E-p 

A = A^+,A^ ^ ( a^_^ b g ) = ^2) 



V = An+lVn = 



E~f} 
ah 



V 




u. 



(C3) 



Let us calculate trM and tr{AP). We find 

trM = tr{A) + tr{B) ^ 2P - 2R (C4) 
tr{M^) = tr{A'^) + tr{B^) + tr{UV) + tr{VU) 

= 2(P + Q)2 + 2(i?2 _2), (C5) 

where 



trM = 2 cos kN. 



(B14) 



The density of states (DOS) D±{E) is calculated for each 
channel ±, respectively: 



dk± = ^dcos~^[x±(£;)/2] 



(B15) 



4-a;l 



-.dE = D±{E)dE. (B16) 



N 



Therefore, the total DOS is given as the sum of D+{E) 
and D_{E): 



D{E) ^D+{E)+D.{E), 



(B17) 



ab 



ah 



ah 



The discriminant D is given by £> = 2tr{AP) - {trM^ + 
8 = 4(5^. As a result, we obtain 

2cos2A: = ^{trM±^/D) (C7) 
= P-R±Q (C8) 

ah ah ah 

Solving the above for E, we can obtain the energy bands. 



20 



.(±) 
^(±) 



^({a + /3)-v± v/(a - 13)^ + + 4(a2 + P) + 2{a + + 4/3?; + 8a&cos 2fc) 
^(^{a + /3)+v + ±^/{a - /3)2 +1-2+ 4(a2 + 62) - 2(a + l3)v - 4/3w + 8a6cos2fc 



(CIO) 
(Cll) 



In Fig. 1(a), the energy band structure for a = b = 1 
when a — I, (3 = is given with varying the interchain 
hopping V. Figure 1(b) shows the cross-section view at 



Ek 




FIG. 18: Energy bands of the decorated ladder model, (a) 
The energy bands as a function of V. (b) The snap shot of 
the energy bands when V = 1. k means the wave vector in 
units of J such that —1.0 < A: < 1.0, V means the 7r-electron 
hopping integral between the inter chain sites, and E means 
the energy in units of = 1. Here we have taken the values 
a = 1.0, /3 = 0, a = 6= 1. 



The extension to the decolated three chains (d = 3) is 
straightforward. Here we give the transfer matrix for a 
simple system of coupled three chains. = A^n$n 

where = (</);^+i, (/)^, 0^, 0^)*. 6 x 6 trans- 

fer matrix with the 2x2 submatrices is given in block 
tridiagonal form: 





* n 





^ n 









^ n 


Bn 



(C12) 



where 



1 







' ^ n 



A„ + i,„ 



u„ 





' ^ n 








For simplicity, we set A„4.i 



n+l,n 



til Cn+1,'. 



a, Cn 



K, and Vn 



Un 



V at site n. 
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Even for general multichain models, some useful for- 
mula exist in order to obtain the eigenvalues of block 
tridiagonal matrices and the determinants of the corre- 
sponding block tridiagonal matrices (68l - [70| . 

Appendix D: Nonlinear Schrodinger Equation 

In this appendix, we derive the discrete nonlinear 
Schrodinger equation for ID electronic system coupled 



Here m„ and f „ are lattice displacements and the internal 
vibration coordinates of the nth unit cell and p„ and P„ 
are their conjugated momentum, respectively. And M is 
the mass of the unit cell, Wo is the oscillation frequency 
of the optical phonon and the dispersion relation of the 
acoustical motion is LUsik) = Cgk, where Cs is the sound 
velocity along the chain. Note that the two kinds of os- 
cillations can be regarded as the dynamics of radial and 
angular coordinates in the polaron models in Sect. 4. 

Then in the Hamiltonian for electrons, both the on- 
site potentials En and the hopping integrals Vm,n depend 



where a and 7„ are the coupling constants. Here Vb is the 
bare amplitude of the hopping term and the term 5Vnn+i 
is a random contribution from the conformational disor- 
der, which we include to describe. Anderson localiza- 
tion takes place when i?„ is a static on-site randomness. 
When En = constant in the system with half-filled con- 
duction band, such as an ID ionic crystal, the SSH term 
generates dimerization in the ground state (the Peierls 
instability) and forms solitons in the excited states. 

However, DNA is considered as a band insulator. Both 
interactions can generate lattice distortions and lead 
to polaron formation when a charge is doped into the 



with lattice oscillations. The relative motions between 
two different base pairs can be represented by an acous- 
tical phonon mode and the vibrational motion inside a 
base pair by optical phonons. These modes represent the 
lattice distortions such as sliding, twisting or bending. 
The Hamiltonian that describes these modes is given by 



(Dl) 



I 

upon these vibrations, in principle. The charge coupling 
to the acoustical phonons is given by the Su-Schrieffer- 
Heeger (SSH) model [13, [EH and the interaction with 
the optical phonons is described via the molecular crys- 
tal model of Holstein [7lj. The SSH model deals clas- 
sically with the lattice degrees of freedom, while elec- 
trons are treated quantum mechanically. Thus the total 
Hamilonian for the electron-phonon interactions in the 
DNA system is given by 



(D2) 



I 

molecule. 

The coupling with (u„+i — u„) usually induces small 
polarons in ionic crystals. The calculations showed that 
a polaron may be built and be robust within a wide range 
of model parameters. The influence of random base se- 
quences was apparently not strong enough to destroy it. 
Thus, polaron drifting may constitute a possible trans- 
port mechanism in DNA oligomers. 

From the total Hamiltonian Htot = Etph + E^ei-ph of 
the system we can derive the equations of motion for the 
variables </>„, u„ and u„ as follows: 



^^2M 2 ° 



+ Un+l - Un + 5Vnn+l){.CiCn+l + Cl^^Cn) + ^(^^n + 7n«n)CtC„, 
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+ Ujl{2Un{t) - Un+l{t) ~ Un-l{t)) - 2aRe{(l>*^+,{t)(f>n{t)(bn *{t)^n~i{t)) (D3) 

= (-En +7„w„)0„(t) +^im„0„(i), (D5) 

m 

I 



where tmn = (Vb + u„i — u„ + (5T/„i„). This means that 
the equihbrium position of each atom in the lattice is 
charged by an amount proportional to the probability 
for the electron to occupy that special atom. The twist 
polaron modify the inter-base electronic coupling, though 
this effect is apparently less strong than the coupling in 
the Holstein model. Accordingly, we also assume that 
the dependence of the hopping integrals on the u„ is so 
week so that it can be ignored. In the following part, we 
deal with only the Holstein polarons. 

Additionally, we assume that time scale of lattice vi- 
brations and electron evolution are as such that vibra- 



tions are slaved by electron probability. It is then possi- 
ble to set time derivative in Ea. (|D4[) equal to zero. Then, 
we obtain 

= W0m(t) + a|(/.„(<)|2</,„(t) + £;„0„(i)(D6) 



where a — By setting t„i„ — —t if |m — n| — 

1 and zero otherwise, we obtain the discrete nonlinear 
Schrodinger equation (DNSE), 



'^-^ = -K^n+iit) + K-iit)) + a|0„(i)|'0n(i) + EMt). (D7) 
at 

I 



The time-independent version is given as 

E(l)n = -t{(j)n+l + (j)n-l) + a\(j)n\^ + E^cjin, (D8) 

, which appears in the Holstein polaron model on the lat- 
tice. The DNSEs (40) and (41) show more various prop- 
erties due to the inhomogeneity of En and the strength 
a for nonlinearity, etc ^73] . 

For example, localization-delocalization transition 
takes place, depending on the coupling strength and the 
initial state. The DNSE has been studied in the context 
of delocalization due to nonlinearity. Indeed, Eq. (jD7p 
describes the ID disordered waveguide lattice, which 
is called the Gross-Pitaevsky (GP) equation on a dis- 
cretized lattice [zl] . When the absolute value of the non- 
linearity parameter a is greater than some critical value 



Qfc, the excitation is self-trapped. It was mathematically 
proved that the DNLS has quasi-periodic self-trapped so- 
lutions called the discrete breathers ^75, 76]. On the other 
hand, it is found that at moderate strength of nonlinear- 
ity the spreading of the wavepacket algebraically grows 
as V(An)2) ~ t''{^Ji ~ 0.2) [t^. 

Appendix E: Coupled Nonlinear Schrodinger 
Equations 



In this appendix, we derive the coupled nonlinear 
Schrodinger equations in order to describe the effect of 
formation of a double strand of DNA on the polarons in 
the HOMO band, following the Holstein's argument [l2| . 
The energy of HOMO band is expressed as, 



N 



i^ff ({x„ , 2/„ }) = Fpft ({x„ , 2/„ }) + ^ (eff (a;„) 1 0;^ 1 2 -f (y„) 1 0f 1 

n=l 

N N 



(El) 



Differentiating the energy with respect to Xp, yp, respec- ordinates Xp, jjp by = 0, — 0; 
tively, we approximately obtain the most contributed co- ^ „ . „ ^ . 

X, = (E2) 

A-B|2 ,,i'j,AiB* , ±B ,A* 



23 



where p = Fh/IquJq, v 
we have assumed u„ — 



= a//otiJo. In the derivation, 

V ~ a\xp ~ = 0, 



dSi 



dxp 



^ = 0, where S^^ = + 0^0;^*) 

I'/'n Pd'/'n P) nieans the frontier orbital density at p-th nu- 
cleotide group in chain A(B) and 5'^'^ means the over- 
lapping integral of the frontier orbitals for electrons in 
the HOMO at p-th nucleotide group between the chains 
A and B. Substituting the expressions into Eq.(34) and 
(35), we obtain the following coupled DNSEs: 



AB iB 



n+1 



1-1/ + 



BAiB 



.ABiA 



E<t>i (E4) 
E^^ (E5) 



where 



AB 
BA 
„.AB 



A\2 



B\2 



pFh\ 

4-pFh\ 



AB 



B\2\ 



(E6) 
(E7) 
(E8) 



By the same argument, similar equations can be obtained 
for the LUMO band case as well, just by replacing the 
orbitals of electrons in HOMO bands with the frontier 
orbitals for holes in the LUMO bands, respectively. In 
the limit w — >■ 0, — )■ 0, it becomes the decoupled DNSE 
in Eg.ljDSp without randomness. 



Appendix F: Lyapunov Exponents and Multichannel 
Conductance 



The definition for energy dependence of the Lyapunov 
exponents is given by, 

7, = hm -l^loga,(Af,(7V)tM,(A^)), (Fl) 

where ai{...) denotes the i-th eigenvalue of the 
matrix Md{N)Hld{N) [s^. As the transfer ma- 
trix Td{N) is symplectic, the eigenvalues of the 
Md{Ny Md{N) have reciprocal symmetry around the 
unity as e'^^ , e'^'^e^'^'' , e"'^^ , where 71 > 72 > 
■ ■ - Id > 0. The d denotes the number of channels; i.e. 
d = 2 in the two-chain model and d = 3 in the three-chain 
model. 

The Lyapunov exponent is related to the DOS p{E) 
as an analogue of the Thouless relation (the generalized 
Thouless relation) ITSI: 



center depend on the parity, the bipartiteness and the 
boundary condition. Since discussions on the details is 
out of scope of this paper, we give simple comments. 
Note that the parity effects appear in the odd number 
chains with the hopping randomness. In the odd number 
chain with the hopping randomness, only one mode at 
_E = is remained as an extended state, i.e. 7^ = 0, 
while other exponents are positive, 7^-1 > • • ■ > 71 > 0. 
The behavior is seen in Fig. 14 in the main text. Then, 
non-localized states with 7 = determine the conduc- 
tance. Although we have ignored the bipartite structure 
in the three-chain models for simplicity, if we introduce 
bipartiteness in the intrachain hopping integral Vn{= Un) 
another delocalized state due to chiral symmetry appears 
at = 0. 

Furthermore, we find that in the thermodynamic limit 
[n — >■ cxd), the largest channel-dependent localization 
length S^d — ^/id determines the exponential decay in the 
Landauer conductance g{n) that is measured in units of 

/h at zero temperature and serves as the localization 
length of the total system of the coupled chains [13, |4l| . 
It is given as 



1 



1 ?.(n) 



1 



exp(- 



2n , 



(F3) 



for n — > 00, where n denotes the system size along the 
chain [3^. Recently, electron transport for molecular 
wires between two metallic electrorodes has been also 
investigated by several techniques. 

If we impose a voltage difference V over the molecular 
structure, we have one end point with energy Ep + eV 
and the other still Ep- Only electrons in that range con- 
tribute to the conductance. Therefore, the I —V charac- 
teristics can be evaluated as 

(F5) 



where f{E) = l/(e^/'^«'^ 4 
tion function and Tij = \tij 



1) means the Fermi distribu- 
P is the squared transmission 
amplitudes between the i-th and j-th channels. At the 
zero-temperature limit, all quantities can be evaluated 
by the transmittion matrix at E ^ Ep. 



d r, 

T.^^{E)- J 



In \p{E )dE 



(F2) 



Accordingly, we can see that the singularity in the largest 
Lyapunov exponent is strongly related to the singularity 
in the DOS. 

Generally speaking, in the quasi-one-dimensional chain 
with the hopping disorder, the singularity in the DOS, 
the localization length and the conductance at the band 



Appendix G: Modified Bernoulli Map 

The correlated binary sequence {Vn} and/or {C„„+i} 
of the hopping inte gra ls can be generated by the modified 
Bernoulli map 

mm-- 



" ^ x„ - 2^^-1(1 - {Xn e h), 



(Gl) 
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where Jo = [0,l/2),/i = [1/2,1). Bq and Bi are 
the bihircation parameters that control correlation in 
the sequence. We set 1 < Bq < -Si < 2 for sim- 
plicity. The asymmetry of the map {Bq ^ B\) cor- 
responds to the asymmetric property in the distribu- 
tion for a real sequence in the double helix DNA, where 
the number of the A-T pairs is not equal to that of 
G-C pairs; they are different from random binary se- 
quences with equal weight. We introduce an indicator 
Rqc for the rate of the G-C pairs in the sequences such 
as Rgc = {Ng + Nc)/{Ng + Nc + Na + Nt), where 
Ng, Nc, Na and Nt denote the numbers of symbols G, 
C, A and T in the base sequence, respectively. 

In the case of Bq = -Bi(= B), depending on the value 
B, the correlation function C(n)(= (KioKio+n)) (^o = 1; 
n is even ) decays following inverse power-law as C(n) 

— 2 — B 

n for large n (3/2 < B < 2). The power spectrum 

2B-3 

becomes S{f) ^ f ^-i for small /. We focus on the 
Gaussian and non-Gaussian stationary regions {1 < B < 
2) that correspond to some real DNA base-pair sequences 
with S{f) - /-"(0.2 <a< 0.8). 

In the ladder models of Sect. 5 , we use for the in- 
terchain hopping integrals at odd sites n the symbolized 
sequence {Vn} that are defined by the following rule: 



In the three-chain models, the interchain hopping in- 
tegral Vn = Un for every site n can be generated in the 
same way tas to the two-chain models. Furthermore, we 
use the successive sequence {Ar„,Ar„+i}, when we make 
a correlated binary sequence {C„„+i} as the hopping in- 
tegral of the middle (nucleotide) chain as follows: 



a 



nn+1 



WaT-AT {X„ € Io,Xn+l e lo) 

Wat-gc (x„ e /o,^n+i e /i) 

Wgc-at {Xn e e h) 

WgC-GC {XnG h,Xn+lG h). 



(G4) 



In the numerical calculation, we assume the following 
rules for simplicity: 



Wgc~gc = Wat-at 1"^ 

Wat-gc = (Wat-at + Wgc-gc)/'^ (G5) 
Wgc-at = Wat-gc- 
As a result, the parameters are Wat and Wat -at- This 
simple rule is based on the binary classification of the 
four DNA nucleotides: adenine and guanine are purines; 
cytosine and thymine are pyrimidines. 



Vn 



Wat = Wta {Xnelo) 

WgC = WcG {Xn&h). 



(G2) 



(Sec Fig. 10(a).) In the mmicrical calculation, Wgc is set 
at a half of Wat (Wgc = Wat/2) for simplicity. Then, 
the artificial binary sequence can be roughly regarded 
as the base-pair sequence as observed in the A— DNA. 
Note that the correlated four-letter sequence {En} such 
as DNA can be also generated by using two independent 
sequences {Xn}, {Yn} by the modified Bernoulli map: 



En = < 



A (X„e Jo,F„€/o) 

T (X„G/o,r„G/i) 

G (x„e/i,y„e/o) 

C (x„e/i,y„e/i). 



(G3) 
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